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ABSTRACT 

We report on multi-wavelength measurements of the accreting black hole Swift J1753.5-0127 in the hard 
state at low luminosity (L ~ 2.7 x 10 36 erg s _1 assuming a distance of d = 3 kpc) in 2014 April. The radio emis¬ 
sion is optically thick synchrotron, presumably from a compact jet. We take advantage of the low extinction 
(E(B-V) = 0.45 from earlier work) and model the near-IR to UV emission with a multi-temperature disk model. 
Assuming a black hole mass of Mbh = 5M 0 and a system inclination of i = 40°, the fits imply an inner radius 
for the disk of R{ n /R g > 212 c/3 (Mbh! 5Mq) _ i , where R g is the gravitational radius of the black hole, and is 
the distance to the source in units of 3 kpc. The outer radius is R oa t/R g = 90.000c/3 {Mbh! 5M 0 ) _1 , which cor¬ 
responds to 6.6 x 10 10 11 z/3 cm, consistent with the expected size of the disk given previous measurements of the 
size of the companion’s Roche lobe. The 0.5-240 keV energy spectrum measured by SwiftIX RT, Suzaku (XIS, 

PIN, and GSO), and NuSTAR is relatively well characterized by an absorbed power-law with a photon index 
of T = 1.722 ± 0.003 (90% confidence error), but a significant improvement is seen when a second continuum 
component is added. Reflection is a possibility, but no iron line is detected, implying a low iron abundance. 

We are able to fit the entire (radio to 240 keV) spectral energy distribution (SED) with a multi-temperature 
disk component, a Comptonization component, and a broken power-law, representing the emission from the 
compact jet. The broken power-law cannot significantly contribute to the soft X-ray emission, and this may be 
related to why Swift J1753.5-0127 is an outlier in the radio/X-ray correlation. The broken power-law (i.e., the 
jet) might dominate above 20 keV, which would constrain the break frequency to be between 2.4 x 10 10 Hz and 
3.6 x 10 12 Hz. Although the fits to the full SED do not include significant thermal emission in the X-ray band, 
previous observations have consistently seen such a component, and we find that there is evidence at the 3.1-cr 
level for a disk-blackbody component with a temperature of kT m = 1 SOl^ eV and an inner radius of 5-14/C,. 

If this component is real, it might imply the presence of an inner optically thick accretion disk in addition to 
the strongly truncated (R[ n > 212 R g ) disk. We also perform X-ray timing analysis, and the power spectrum is 
dominated by a Lorentzian component with u max = 0.110 ±0.003 Hz and // rnax = 0.16 ±0.04 Hz as measured by 
XIS and XRT, respectively. 

Subject headings: accretion, accretion disks — black hole physics — stars: individual (Swift J1753.5-0127) 

— X-rays: stars — X-rays: general 
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1. INTRODUCTION 

Most accreting stellar-mass black holes in binary sys¬ 
tems exhibit large changes in luminosity over time, ranging 
from a substantial fraction of the Eddington limit (L K dd) to 


~ 10“ 8 or ~ 1 0 _y Z-Kdd- In addition to changes in luminos¬ 
ity, these systems show other observational changes, includ¬ 
ing transitions between distin ct sp ectral states that are similar 
from system-to-system (e.g., McClintock & Remillardl 120061 : 
iBelloniluOlOt) . The thermal dominant (or soft) state has a 
strong thermal component from an optically thick accretion 
disk in the X-ray spectrum. In the hard state on the other hand, 
this component contributes a lower fraction of the flux in the 
X-ray band. The drop in flux is partly due to a decrease in the 
temperature of the component dKalemci et al.ll2004l) . moving 
its peak into the ultraviolet where it is difficult to measure due 
to interstellar absorption. While the soft thermal X-ray emis¬ 
sion weakens, there is a strong increase in the hard X-rays, 
and the X-ray spectrum in the hard state is dominated by a 
power-la w, which often has an exponential cutoff abo ve 50- 
100 keV dGrove et alJll998tlGilfanov & MerlonilHoH) . 

Accreting black holes also emit in the radio band when 
they are in the hard state, and this is due to a powerful com¬ 
pact jet dCorbel et alJl200(i lFendetll2001l) . At radio frequen¬ 
cies, the spectrum is dominated by a partially self-absorbed 
synchrotron component that has a flat or rising spectrum 
(F v oc ;/\ where a )> 0). The jet spectrum changes slope 
above the break frequency, i^ak, becoming steeper because 
the frequency is sufficiently high that the entire jet is optically 
thin. In some cases, the measurement of rv ea k has been con 


strained to be in the infrared (IR) to optic al (Corbel & Fender 
2002^ Gandhi et alJ 1201 lb iRahoui et al] 1201 It iRussell et al 


201 3al bl 120141) . but its measurement can be complicated be 


cause of the other emission components (e.g., from the ac¬ 
cretion disk or the optical companion) and also because the 
jet spectrum is likely significantly more complicated than a 
simple broken pow er-law (iMarkoff, Nowak & Wilmsl 12005b 
iMigliari et alJl2007l) . 

In the hard state, it is clear that there is a strong con¬ 
nection between the X-ray and ra dio emission. The fluxes 
in the two bands are correlated (tCorbel et alJ l2000l 120031: 

Corbel, Koerding & Kaaret 

120141) . and while early 


Gallo, Fender & Pooled 

200 j iCorbei et alJl2013t 


2003: 


Gallo et al. 


studies suggested that all black hole sources might lie on 
the same correlation line, obse rvations of more systems have 
show n that this is not the case dJonker et al.l 2010 l : ICoriat et alJ 
1201 ll) . A current topic of debate is whether all sources 
lie on two correlation lines, one track for standard sources 
and one f or outliers, or if there is a continuum of differ- 
ent tracks (ICoriat et al.ll201 it ICorbel et alJl2013t iGallo et aP 
120141) . Another topic is how much, if any, of the X-ray 
emission originates in the jet. While the most typical hard 
state spectrum with an exponential cutoff is well described 
by thermal Comptonization, and it has been argued that it 
is unlikely that this_emission is due to synchrotron emis¬ 
sion from a jet (Zdziarski et al. 2003), some black hole spec¬ 
tra ap pear to have multiple high-energy continuum compo¬ 
nents dJoinet et al.l2007HRodriguez et alJl2008l:lBouchet et all 


2009 : iDroulans et al.l l2Q10i: Russell et alJ 12013) 7 and a jet 
origin is not ruled-out. In fact, Cygnus X-l often shows 
two high-energy com ponents in the hard st a te, including 
an MeV component dMcC onnell et alJ 12 000: Ra houi et alJ 
1201 ll IZdziarskkF ubihski & Sikorall2012l) . and the detection 
of strong polarization at >400 keV favors a synchrotron ori¬ 


gin dFaurent et al.ll20ll1 : lJourdain et al.ll2012l) . 

The fact that a jet is present in the hard state and that there 
is some connection between the disk and the jet leads to the 
question of what we know about the disk properties. The main 
question regarding the optically thick disk concerns the loca¬ 
tion of the inner radius (R m ). One idea is that the black hole 
states are essentially determined by the mass accretion rate 
and R m (lEsin. McClintock & Naravanll997l) . with sources en¬ 
tering the hard state because of an increase in Rj n . How¬ 
ever, X-ray observations of sources in the bright hard state 
seem to contradict this since relativistic ally smeared reflec¬ 
tion components are seen from some systems that imply that 
the dis k remains close to the innermost stable circular orbit 
(ISCO:lBlum et aljl2003 IReis et alJl20lH IMiller et akllloH 

Fabia n et al.l 120121) . In addition, ther mal component mod- 

eling has led to similar conclusions (IReis, Fabian & Milled 
12010 ). While photon pile-up in CCD spectra has_sparked 
some debate about iron line results (IMiller et al.l 120061 : 
iDone & Diaz Trigol2010HMiller et al.l2Q10l) . more recent ob¬ 
servations with the Nuclear Spectroscopic Telescope Array 
(NuSTAR) confirm strongly broadened and skewed iron lines 
in the bright hard state for GRS 1915+105 (IMiller et al | 20 1 3|). 
GRS 1739-278 (IMilleretal.ll 20 Tl). GX 33 9-4 dFurst etal] 
120151) . and Cvgnus X-l dParkeretal] 120151) . For the case 
of GRS 1739-278, the luminos ity is ~5 % Lf,im, and the 
inferred inner radius is < 1 2 R g (IMiller et alJ 120151) . where 
R g = GM BH /c 2 and G is the gravitational constant, Mbh is 
the black hole mass, and c is the speed of light. Signifi¬ 
cantly truncated disks have been reported for the hard state 
at intermediate a nd low luminosities using reflection com¬ 
ponent modeling dTomsick et alJ 120091 : IShidatsu et alJ 1201 ll : 
iPlant et aD 120151) and also by mod eling the thermal compo¬ 
nent from the optically th ick disk (Gierlin ski, Done & Page 
[20081: iCabanac et al.ll2009l) . 

To investigate questions related to the accretion geometry 
and the relationship between the disk and the jet, we per¬ 
formed multi-wavelength observations of the accreting black 
hole Swift J1753.5-0127 in the ha rd state. This syste m was 
first discovered in outburst in 2005 dPalmer et al.l2005l) . and it 
is very unusual in that it has been bright in X-rays for almost a 
decade. The optical light curve shows a 3.2 hour modulation, 
which has been interpreted as a superhump period (a modula¬ 
tion due to tidal stresses on a precessing, elliptical accretion 
disk), suggesting that the orbital period is somewhat smaller 
than this (Zurita et al. 2008). From radial velocity measure¬ 
ments, iN^^^^veTdJ^Ol^) find a 2.85 hour signal, which 
is likely the true orbital period. Thus, Swift J1753.5-0127 has 
one of the shortest orbital periods of any known black hole bi¬ 
nary. Although the mass of the black hole in Swift J 1753.5— 
0127 is still debated since there has not been an opportunity to 
obtain a radial velocity measur ement for the compani on star 
with the system in quiescence, iNeustroev et alJ d2014i) argue 
that the mass is relatively low, Mbh < 5M 0 , and we adopt a 
black hole mass of 5 M@ for calculations in this paper. 

Swift J1753.5-0127 is also unusual in that it has a low 
level of extinction, due in part to it being somewhat out 
of the plane with Galactic coordinates of / = 24.9° and b = 
+12.2°. iFroning et alJ (120141) obtained UV measurements 
showing that E(B-V) = 0.45, and we confirm this value in 
a companion paper (Rahoui et al., submitted to ApJ). It is 
not entirely clear whether the system is relatively nearby or 
in the Galactic halo as there is a large range of possible 
distances, d = l-10kpc dCadolle Bel et al.1120071: iZurita et alJ 
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l2008t iFroning et al.ll2014l) . iFroning et al.l (120141) provide evi¬ 
dence that the UV emission from Swift J1753.5-0127 in the 
hard state comes from an accretion disk, and they calculate 
distance upper limits that depend on Mbh, assuming that the 
mass accretion rate is less than 5% LecM- For AT BH = 5M 0 , 
the upper limit is 2.8-3.7kpc, depending on the inclination 
of the system, and we use a fiducial distance of 3 kpc for the 
calculations in this paper. 

Swift J1753.5-0127 has been extensively observed in the 
radio band, and is one of the clearest examples of a source that 
is an outlier in the radio /X-ray correlation plot (ISoleri et al.1 
120101: ICorbel et alJl2013h . The location on the plot depends 
on the assumed distance, and the previous work has assumed 
a source distance of 8 kpc. While we are adopting a signif¬ 
icantly smaller distance, ISoleri et akl (1201 Ol) considered how 
distance affects the the radio underluminosity, which is a 
measure of how far a source is from the standard correla¬ 
tion. ISoleri et al. (2010 ) show that a smaller distance moves 
the source farther from the standard correlation (see also 
iJonker et al.ll2004 ). Thus, the fact that recent work suggests 
that Swift J1753.5-0127 is closer than early estimates only 
strengthens the conclusion that the source is an outlier. 

For this work, we have carried out a large campaign to ob¬ 
serve Swift J1753.5-0127 in the hard state with radio, near- 
IR, optical, UV, and X-ray observations as described in Sec¬ 
tion 2. In the X-ray, data were obtained with NuSTAR, Suzaku, 
and Swift/X RT. The observations occurred when the flux level 
was close to the minimum brightness this source has had in 
the ~ 10 years since its discovery (see Figure Q}. The low 
flux level (and presumably mass accretion rate) may cause 
changes in the properties of the accretion disk or jet compared 
to previous observations at higher flux levels. In Section 3, we 
perform spectral analysis for the different energy ranges (ra¬ 
dio, near-IR to UV, and X-ray) separately and then also as a 
combined radio to X-ray Spectral Energy Distribution (SED). 
We also produce an X-ray power spectrum for timing anal¬ 
ysis. We discuss the results in Section 4, and then provide 
conclusions in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations that we obtained in 2014 April are listed 
in Table Q] and more details about the observation times are 
shown in Figure[2] The X-ray flux was rising very slowly dur¬ 
ing the observation, and this is seen especially clearly in the 
Suzaku/XIS light curve (Figure EJ. We provide more details 
about the observatories used and how the data were processed 
in the following. 


2 .1. Radio 

We observed Swift J1753.5-0127 with the Karl G. Jan¬ 
sky Very Large Array (VLA) on 2014 April 5 (MJD 56752) 
from 11:00-13:00 UT with the array in its most-extended A- 
configuration. We split the observing time between the 4— 
8 and 18-26-GHz observing bands. In the lower 4-8 GHz 
band, we split the available bandwidth into two 1024-MHz 
basebands, centered at 5.25 and 7.45 GHz. Each baseband 
was split into eight 128-MHz spectral windows, each of which 
comprised sixty-four 2-MHz channels. The higher-frequency 
18-26 GHz band was fully covered by four 2048-MHz base¬ 
bands, each comprising sixteen 128-MHz sub-bands made 
up of sixty-four 2-MHz channels. After accounting for cal¬ 
ibration overheads, the total on-source integration times for 
Swift J1753.5-0127 were 25.3 minutes in the 4-8-GHz band 
and 29.1 minutes in the 18-26-GHz band. 



MJD-55000 (days) 


Fig. 1.— MAXI (a) and Swift !BAT (b) light curves in the 2-20 keV 
and 15-50 keV bands for Swift J1753.5-0127 between mid-2009 
and mid-2014. The vertical lines mark the start and stop times of 
the observations used in this work (2014 April 2-8). 

The data were reduced using version 4.2.0 of the Co mmo n 
Astronomy Software Application (CASA; iMcMullin et alJ 
12001 , We applied a priori calibration to account for up¬ 
dated antenna positions and gain variations with changing 
elevation or correlator configuration, and corrected the 18- 
26-GHz data for opacity effects. We edited out any data af¬ 
fected by antenna shadowing before Hanning smoothing the 
data and removing any radio frequency interference. At all 
frequencies we used 3C286 to calibrate the instrumental fre¬ 
quency response, and to set the amplitude scale according to 
the default Perley-Butler 2010 coefficients implemented in the 
CASA task SETJY. We used J1743-0350 as a secondary cal¬ 
ibrator to determine the time-varying complex gains arising 
from both atmospheric and instrumental effects. 

The calibrated data on Swift J1753.5-0127 were averaged 
by a factor of four in frequency to reduce the raw data volume, 
and then imaged using Briggs weighting with a robust param¬ 
eter of 1 to achieve the best compromise between sensitivity 
and sidelobe suppression. When imaging, we used the multi¬ 
frequency synthesis algorithm as implemented in CASA’s 
CLEAN task, choosing two Taylor terms to account for the 
frequency dependence of source brightness. The source was 
clearly detected in all frequency bands, with an inverted spec¬ 
trum (a > 0, where the flux density varies as a function 
of frequency v as S v oc v a ). To better constrain the radio 
spectrum, we split each frequency band into four frequency 
bins (of width 1024 MHz at 4-8 GHz, and 2048 MHz at 18- 
26 GHz where the intrinsic sensitivity per unit frequency is 
lower). We measured the source brightness in each frequency 
bin by fitting an elliptical Gaussian to the brightness distribu¬ 
tion in the image plane. Swift J1753.5-0127 appeared unre¬ 
solved at all frequencies. 

Swift J1753.5-0127 was also observed with the Mullard 
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Radio Astronomy Observatory^ Arcminute Microkelvin Im¬ 
ager (AMI) Large Array (IZwart et al.l 120081 ) on two occa¬ 
sions during the coordinated multi-wavelength campaign be¬ 
tween 2014 April 4-5. These ^4-hour observations were 
taken at the times given in Table Q] with a central frequency 
of 15.4 GHz. The AMI Large Array consists of eight 13-m 
dishes, with the full frequency band of 12-17.9 GHz being di¬ 
vided into eight 0.75 GHz bandwidth channels. Channels 1-2 
and 8 were ignored due to lower response in those frequency 
ranges and the expected high level of interference from satel¬ 
lites due to the low elevation of the source. The primary beam 
FWHM is ~6 arcmin at 16 GHz. 

The AMI data were reduced using the semi-automated 
pipeline procedure described in IStalev et~ahl (120131 ). which 
uses the AMI software tool REDUCE to automatically flag 
for interference, shadowing and hardware errors, calibrate 
the gain, and synthesize the frequencychannels to produce 
visibility data in uv-FITS format (see IStalev et al.ll2.013l for 
more details). However, the low elevation and the radio¬ 
quiet nature of the source resulted in high noise levels in 
the reduced images, and thus the two observations were con¬ 
catenated to maximize the signal-to-noise. The concatenated 
dataset was then imaged in CASA, where the CLEAN task 
was used to produce the combined frequency image, and the 
flux density was measured by fitting a Gaussian model to 
the source in the radio map using the MIRIAD task IMFIT. 
The error on the concatenated flux density was calculated as 

a = J(0.05S U ) 2 + Og t + of ms following lAMI Consortium et al.1 

rt 20 l 2 h , with a 5% absolute calibration error added to the fit¬ 
ting error crpt calculated in IMFIT. The source concatenated 
flux density was measured at 290 ± 50 pJy. 

2 .2. Ground-Based Optical and Near-IR 

Kanata is a 1.5 m telescope at the Higashi-Hiroshima 
Observatory. Photometric observations were performed for 
this study on three nights (MJD 56749, 56751, and 56754) 
with t he B , V, 7, and K x bands using the HO NIR instru¬ 
ment (Sakimot o et al.ll2012t lAkitava et al.1120141) attached to 
Kanata. The individual frame exposure times were 75, 136, 
120, and 60s in B , V, /, K s -bands, respectively. The data 
reduction was performed in the standard manner: the bias 
and dark images were subtracted from all images, and then 
the images were flat-fielded. The magnitudes of the ob¬ 
ject and comparison stars were measured using PSF pho¬ 
tometry. For the B , V, and /-band photometry, we used 
the comparison star located at R.A.= 17 h 53 m 25. s 275, Decl.=- 
01°27 , 30."05 (J2000.0), w hich has magnitudes of B = 17.62 
V = 1 6.66, and / = 14.468 (IZurita et al]l2008HSkrutskie et all 
2006). For the A)-band photometry, we used the comparison 
star at R.A.= 17 h 53 m 25 .^853, Decl=-01°26 / 1 7/'00 (J2000.0), 
for which K s = 11.132 (ISkrutskie et alj|2006h . 

We also conducted optical g ', r', i', and z' band monitor¬ 
ing observations with the Lulin 41 cm Super-Light Telescope 
(SLT), which is located in Taiwan, on three nights in 2014 
April (see Table Q}. Photometric images with 180 s exposures 
were obtained using the U42 CCD camera. We performed the 
dark-subtraction and flat-fielding correction using the appro¬ 
priate calibration data with the IRAF package. Photometric 
calibrations were made with the Pan-STARRSl 37 t catalog s 
dMagnier et al.ll2013l : ISchlaflv et al .1120 1 2 t iTonrv et ahll2012l ) . 

The DAOPHOT package was used to perform the aperture 
photometry of the multi-band images. 



MJD-55000 (days) 

Fig. 2.— (a) The filled circles show the MAXI light curve over the 
time of the observations (2014 April 2-8). The times of the AMI, 
VLA, 5w(/f/XRT, Suzaku, and NuSTAR observations are indicated. 
Also, the 1-12 keV light curve for Suzaku/X ISO is shown (the actual 
count rate divided by 70). (b) The Svii/f/BAT light curve over the 
time of the observations, (c) The optical and near-IR magnitudes 
measured at the Kanata and SLT telescopes. 

2.3. Swift 

The Swift satellite dGehrels et all [2004) includes 
two pointed instruments, the X-ray Telescope (XRT; 
iBurrows et alJ 120 051) a nd the Ultra-Violet/Optical Telescope 
(UVOT: [Roming et al.) 120051) . and we used data from both 
instruments from ObsID 00080730001 in this work. We 
performed the XRT data reduction using HEASOFT v6.15.1 
and the 2013 March version of the XRT calibration data 
base (CALDB), and made event lists using xrtpipeline. 
The XRT instrument was in Windowed Timing mode to 
avoid photon pile-up. For spectral analysis, we extracted 
photons from within 47" of the Swift J1753.5-0127 position, 
and made a background spectrum from a region away from 
the source. We measured an XRT source count rate of 
7.7 c/s in the 0.5-10keV band during the 2.4ks observation. 
We used the appropriate response file from the CALDB 
(swxwt0to2s6_20010101v015.rmf) and produced a new an¬ 
cillary response file using xrtmkarf and the exposure map 
generated by xrtpipeline. We binned the 0.5-10keV 
spectra so that each bin has a signal-to-noise of 10 . 

For UVOT, we obtained photometry in six filters (v, b, u, 
uvw\, uvm2, and uvw2) during the observation. For each fil¬ 
ter, we produced an image using uvotimsum and made a 
source region with a radius of 5" and a background region 
from a source-free region. Then, we used uvotsource to 
perform the photometry and calculate the magnitude and flux 
of Swift J1753.5-0127 for each filter. 

Swift also includes the wide field of view (FoV) Burst Alert 
Telescope (BAT), and we use data from BAT to study the 
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long-term 15-50 keV flux (see Figures Q] and [2j. 

2.4. NuSTAR 

The Nu clea r Spectroscopic Telescope Array ( NuSTAR ; 
lHarrison et al.l 12013!) consists of two co-aligned X-ray tele¬ 
scopes, FPMA and FPMB, sensitive between 3-79 keV. To 
reduce the data, we used nupipeline v.1.3.1 as distributed 
with HEASOFT 6.15.1. During our analysis, an updated ver¬ 
sion became available, but we carefully checked that it does 
not influence our results. We extracted the source spectrum 
from a circular region with 90" radius centered on the J2000 
coordinates. Due to the triggered readout of the detectors, 
pile-up is not a concern for NuSTAR. The background was ex¬ 
tracted from a circular region with a 170" radius at the other 
end of the FoV. Small systematic changes of the background 
over the FoV can be neglected, as Swift J1753.5-0127 is a 
factor of six brighter than the background, even at 70keV. 
The spectrum includes data from two NuSTAR ObsIDs. We 
reduced both ObsIDs separately and added the resulting spec¬ 
tra and response files using addascaspec. The resulting 
total exposure time is given in Table [I] 

2.5. Suzaku 

For Suzaku , we used data from the X-ray Imaging Spec¬ 
trometers f XISs: [Kovama et aTI 12007) and from the Hard X- 
ray Detector (HXD; Takahash i et akl 20071) PIN diode detec¬ 
tor, and the HXD gadolinium silicate crystal detector (GSO). 
The XIS has three CCD detectors (XIS0, XIS1, and XIS3) 
that operate in the 0.4-12 keV bandpass. We produced event 
lists for each detector using aepipeline and merged the 
event lists taken in the 3x3 and 5x5 CCD editing modes. 
We ran aeattcor2 and xiscoord on each of the merged 
event files to update the attitude correction because this is 
important for the pile-up estimate, which we calculated us¬ 
ing pileest. We extracted source spectra using a 4'-radius 
circle with the inner 22" removed due to pile-up at a level 
of >4% in the core of the point spread function (PSF). We 
extracted the background from a rectangular region near the 
edge of the active area of the detector. The XIS detectors 
were in 1/4 window mode for the observation, and part of the 
source region falls off of the active region of the detector. We 
accounted for this when determining the background scaling. 
We used xisrmfgen and xissimarfgen to produce re¬ 
sponse matrices, and we combined the XIS0 and XIS3 spectra 
(the two front-illuminated CCD detectors) into a single file. 

For HXD, we analyzed both PIN and GSO data us¬ 
ing the Perl scripts hxdpinxbpi and hxdgsoxbpi, re¬ 
spectively, after screening with the standard selection crite¬ 
ria. These scripts produce deadtime-corrected source and 
background spectra automatically. The non X-ray back¬ 
ground model was taken from the FTP siteQ. and cos¬ 
mic X-ray background (CXB) was also subtracted based on 
previous High Energy Astronomy Observatory (HEAO) ob¬ 
servations dGruber et all 11999 4 for PIN. As an energy re¬ 
sponse, we used ae_hxd_pinxinomell_20110601.rspfor PIN 
and ae_hxd_gsoxinom_20100524.rsp with an additional cor¬ 
rection file (ae_hxd_gsoxinom_crab_20100526.arf) for GSO. 
The background count rate is significantly higher than the 
source rate for GSO, but we still clearly detect Swift J 1753.5— 
0127 at a rate of 0.740 ± 0.026 c/s. 

24 ftp://legacy.gsfc.nasa.gov/suzaku/data/background/pinnxb_ver2.2_tuned/ 
and ftp://legacy.gsfc.nasa.gOv/suzaku/data/background/gsonxb_ver2.6/ 



Fig. 3.— (a) Power-law fit to the radio spectrum with T = 0.71 ± 
0.05 (a = 0.29 ± 0.05). The measurements are from VLA (8 points) 
and AMI (1 point at 15.4 GHz), (b) Data-to-model ratio. 


3. RESULTS 

3.1. Energy Spectrum 

We performed all of the spectral fits using the XSPEC 
vl2.8.2 software. For the X-ray spectra, we used instru¬ 
ment response files produced using the HEASOFT software. 
For the radio, ground-based optical and near-IR, and UVOT, 
we determined the flux for each data point and then used 
f lx2xsp to produce spectral files and unitary response ma¬ 
trices that can be read into XSPEC. All spectral fits are per¬ 
formed by minimizing the \ 2 statistic. 

3.1.1. Radio Spectrum 

We fitted the radio points with a power-law model (see Fig¬ 
ure 0. and this provides an acceptable fit with a reduced-^ 2 
(% 2 ) of 0.41 for 7 degrees of freedom (dof). The power- 
law photon index is T = 0.71 ±0.05 (90% confidence er¬ 
rors are given here and throughout the paper unless other¬ 
wise indicated), and this corresponds to a spectral index of 
a = 1 — r = 0.29 ± 0.05 (as mentioned above, a is defined ac¬ 
cording to S v oc v a , where is the flux density). We used 
the XSPEC model pegpwrlw, allowing for the power-law 
normalization to be defined as the flux density at 10 GHz, and 
we obtain a measurement of 256 ± 8 //.ly at this frequency. 

3.1.2. Near-IR to UVSpectrum 

The times of the data taken for the near-IR to UV part of 
the spectrum from Kanata, SLT, and Swift/ UVOT are shown 
in Figure[2] The ground-based (Kanata and SLT) observations 
were taken in five epochs over six nights (see Tables [2] and [3] 
for the exact times of the exposures). As the source is vari¬ 
able from night-to-night and also on shorter time scales, we 
used measurements as close to each other in time as possible, 
while keeping the maximum wavelength coverage. The Swift 
observation occurred between epochs 3 and 4, and we used 
the points from epoch 3 because the Kanata V- and /-band 
measurements occurred on the same night. We also used the 
//-band measurement from epoch 1 because the statistical er¬ 
ror bar is large enough to account for source variability. We 
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did not include the B-band measurement because UVOT cov¬ 
ered the same frequency, and the UVOT measurement was 
closer in time to the other observations. For each SLT band, 
several epoch 3 measurements were made, and for the SED, 
we used the average value. We estimated the uncertainty on 
these points by calculating the standard deviation of the mea¬ 
surements. 

We fitted the near-IR to UV spectrum with a power- 
law model with extinction. The XSPEC extinction model, 
redden, is based on the iCardelli, Clayton & Mathisl (119891) 
relationship. The fit is poor (xt = 10.8 for 10 dof), strongly 
over-predicting the X s -band point. The photon index for the 
power-law is T = 0.2 ± 0.2, but we suspect that this is not 
physically meaningful. A somewhat better fit (although still 
far from being formally acceptable) is obtained by replacing 
the power-law with a blackbody (specifically bbodyrad), 
and this model is a much better match to the spectral slope in 
the near-IR. If E(B-V) is left as a free parameter, xt =8.0 for 
10 dof, and we find E(B- V) = 0.54 ± 0.06, kT = 1.9 ± 0.3 eV, 
and a normalization of Rkm/^io = (1.0±0.1) x 10 12 , where 
Rkm is the size of the blackbody in units of kilometers, and 
d io is the distance to the source in units of lOkpc. Fixing 
E(B-V) to 0.45 gives xt = 7.8 for 11 dof (essentially the 
same quality as the fit with the extinction parameter free), 
kT= 1.51 ±0.05 eV, and R km /r/ 10 = (1.2 ±0.1) x 10 12 . 

The fit is worse with a multi-temperature disk-blackbody 
diskbb model (xt = 13.4 for 11 dof); however, a signifi¬ 
cant improvement is obtained if the outer edge of the disk is 
left as a free parameter. We implemented this by using the 
diskir model (fGierliriski. Done & Pageif2008 u 12009}) . We 
turned off the thermalization in the outer disk (/ out = 0), and 
we set the Compton fraction (L c /Ld) to zero. This model 
gives xt = 3-6 for 10 dof, and the near-IR to UV spectrum 
is shown fitted with this model in Figure [4] For the parame¬ 
ters, we obtain kTi„ = 5± 2 eV for the temperature of the inner 
disk and a value of 1.29^] ^ for log r out , where r out = Rout/Rin, 
and R- m and B out are, respectively, the inner and outer radii 
of the optically thick accretion disk. The diskir normal¬ 
ization, which has the same meaning as the diskbb nor¬ 
malization (Vdiskbb = (Rin,kmMo) 2 /cos/, where R in , km is the 
inner radius in units of kilometers, d io is the distance to the 
source in units of 10 kpc, and i is the inclination of the disk) is 
-Vdiskbb = (9^) x 10 9 . Here, we simply note that this implies 
a very large inner disk radius. We consider the implications 
below in detail after using the same model as a component in 
fitting the full SED. 

None of the fits described above are formally acceptable, 
and there are a few possible reasons for this. Of course, the 
first possibility is that the spectrum requires a more com¬ 
plex model than those we have tried. Secondly, it is known 
that t here is significant variability in this part of the spec¬ 
trum (Zuri ta et al.l [2008: Neu stroev et al.1 120141) . and this is 
also seen in Figure [2] Finally, the largest residuals (see Fig¬ 
ure 01 are in the UV where the extinction changes rapidly. 
Uncertainties in the extinction law and the calibration of the 
broad UVOT photometric bins could also lead to the large 
residuals in this part of the spectrum. 

3.1.3. X-ray Spectrum 

We performed a simultaneous fit to the spectra from all the 
X-ray instruments with an absorbed power-law model, allow¬ 
ing for different overall normalizations between instruments. 
To account for absorption, we used the tbabs model with 



Fig. 4.— (a) Fit to the Kanata, SLT. and UVOT spectra with a multi- 
temperature disk model with outer radius as a free parameter. The 
points are not dereddened, and the model assumes E(B-V) = 0.45. 
(b) Data-to-model ratio. 

IWilms, Allen & McCravl d2000i) abundances and lVerner et all 
(119961) cross-sections. As shown in Table 0 the column den¬ 
sity is Nu = (2.01 ±0.05) x 10 21 cm -2 , the power-law photon 
index is T = 1.722 ± 0.003, and this simple model provides a 
surprisingly good fit with xt = 1-40 for 2143 dof. The resid¬ 
uals (see the data-to-model ratio in Figure [ 05 ) do not show 
any evidence for an iron emission line as might be expected 
if there was a strong reflection component. For a narrow line 
in the 6.4-7.1 keV range, the 90% confidence upper limit on 
the equivalent width is <5 eV, and for a line with a width of 
0.5 keV, the upper limit on the equivalent width is <6 eV. The 
Suzaku/ GSO shows a different slope above «80keV, and we 
added an exponential cutoff using the highecut model. A 
cutoff with E c ut = 66^0 keV and £f 0 id = 218 ^'q 1 keV provides 
a large improvement in the fit to the GSO data, but the overall 
xt only improves to 1.39 for 2141 dof. 

Previous work fitting X-ray spectra of Swift J 1753.5— 
0127 has often shown evidence for a thermal disk- 
blackbody component with an inner disk tem¬ 
perature of kT m = 0.1-0.4 keV when the source 
is _ in the hard state (fMiller. Homan & Miniuttil 


2006; Hiemstra et al. 

2009; IChiang et ah] 2010; 

Reynolds et al. 2010 

Cassatella, Uttlev & Maccarone 

20121 Kolehmainen, Done & Diaz Trigo 2014j). Thus, we 


added a diskbb model to the power-law with an exponential 
cutoff, and the xt improves to 1.29 for 2139 dof (see Table’0. 
While this represents a significant improvement (an F-test 
indicates that the significance of the additional component is 
in excess of 12-er), the temperature is much higher and the 
normalization is much lower than previously seen. Our value 
is Vdiskbb = 3.8^| j compared to values of 1000 reported by 
Reynolds et all (120101) and iCassatella, Uttlev & Maccaronei 
(120121) . A value of Vdiskbb = 3.8 would imply an unphysically 
small inner radius. The equation for the inner radius in terms 
of the gravitational radius is 

R jR g = (0.676 d w f 2 VN^)/mBH/M Q )V^i) , (1) 

where / is the spectral hardening factor (IShimura & Takaharal 
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Fig. 5.— (a) Absorbed power-law fit to the Swift J1753.5-0127 
X-ray spectrum, including Swift/XKT (black), Suzaku/X IS03 (cyan), 
Suzaku/X IS1 (purple), NuSTARI FPMA (blue), NuSTAR/FPMB (red), 
SuzakufPIN (yellow), and Suzaku /GSO (orange), (b) Data-to-model 
ratio. 

119951) . For a distance of 3 kpc, = 5, /= 1.7, which is 

a typical value (Shimura & Takahar al 19951) . and i = 40° based 
on the estimate of lNeustroev et al.1 (120141) . we find R m /R g = 
0.26, which puts the inner radius inside the event horizon. 

Figure [5] shows that there is a small deviation from the 
power-law in the hard X-ray band with the residuals increas¬ 
ing above lOkeV and peaking near 25keV. Although there 
is no iron line, this could still be evidence for a weak re¬ 
flection component or an additional continuum parameter. 
Adding a reflection component to the power-law using the 
reflionx model dRoss & Fabianl 120051 ) provides a signif¬ 
icant improvement in the fit to Xv = 1-27 for 2138 dof. The 
reflection covering fraction (determined by calculating the ra¬ 
tio of the 0.001-1000 keV unabsorbed flux in the reflection 
component to the 0.1-1000 keV unabsorbed flux in the direct 
component) of 0/2 tt = 0.2 and the ionization (parameterized 
by £ = L/nR 2 , where L is the luminosity of ionizing radia¬ 
tion, n is the electron number density, and R is the distance 
between the source of radiation and the reflecting material) 
of £ < 5.3ergcms _1 (see Table |4]) would both be reasonable 
for a cool and truncated disk (although we note that low cov¬ 
ering fractions can also be explained by beaming emission 
away from the disk: iBeloborodovl 19991) . The iron abundance 
of 0.28 ± 0.08 times solar is low but perhaps not unreason¬ 
ably so. Adding a diskbb in addition to reflionx only 
provides a small improvement to the fit (to xt = 1 -26 for 2136 
dof), and Adi sk bb is even smaller than the previous value. How¬ 
ever, it is notable that adding the diskbb component causes 
the iron abundance to change to 0.47!j} 2 j times solar. 

3.1.4. XRT Spectrum and the Possibility of a Thermal Component 

To investigate further on the question of why we do not 
see a physically reasonable diskbb component while many 
previous studies of Swift J1753.5-0127 in the hard state did, 
we fit the X-ray spectra individually. Despite the short ex¬ 
posure time, the Swift/XRT spectrum provides the best infor¬ 
mation on this because it extends down to 0.5 keV without 


strong instrumental features (we note that Suzaku/XIS also 
has sensitivity down at this energy, but the residuals indicate 
the presence of instrumental features). A fit to the XRT spec¬ 
trum with an absorbed power-law model gives Ah = (2.2 ± 
0.2) x 10 21 cm -2 , r = 1.65 ±0.03, and X l= 1-27 for 131 dof. 
Adding a diskbb provides a significant improvement (to 
X 2 = 1.17 for 129 dof), and an F-test indicates a significance 
of 99.8% (3.1-cr) for the diskbb component. The parame¬ 
ters for this fit are Ah = (5 ± 1) x 10 21 cm -2 , F = 1.76 ± 0.06, 
kT m = 130! 2 °eV, and A dlskbb = (1.6!™) x 10 5 . 

Although the column density is not known precisely, it is 
clear that it is lower than «6 x 10 21 cm -2 . The extinction 
value that we use in this paper (E(B-V) = 0.45) corresponds 
to Ah = 3.1 x 10 21 cm -2 based on the relationship derived 
in iGiiver & Ozell d2009t) . Fixing the column density to this 
value and fitting the XRT spectrum with a model consisting 
of a diskbb and a power-law gives thermal parameter val¬ 
ues of kT m = 150^20 eV and ^diskbb = (l-l-as) x 10 4 . Thus, if 
we only had the SwiftIX RT data, we would likely conclude 
that there is a physically reasonable thermal component. The 
kT m = 150eV diskbb component that may be present in the 
Swift/XRT spectrum falls rapidly going to energies below soft 
X-rays and cannot explain the near-IR to UV emission that we 
see. Thus, even if it is real (and it may not be because it does 
not appear to be present when fitting all the available data), 
it is not one of the dominant components in the overall SED, 
and we do not include it in the following as we build a model 
for fitting the full SED. 

3.1.5. Near-IR, Optical, UV, and X-ray Spectrum 

Before fitting the full SED, we fit the near-IR to X-ray 
spectrum in order to determine if it can be fit in a phys¬ 
ically self-consistent manner. As we found that the near- 
IR to UV spectrum requires a thermal model with the outer 
disk radius as a parameter, we start by fitting the spectrum 
with a diskir model. While the fits above used a Comp¬ 
ton fraction of zero (no Comptonization component), here 
we allow L c /L r i to be a free parameter, so that the model in¬ 
cludes Comptonization by a thermal distribution of electrons 
with a temperature of kT e , causing the model to extend into 
the X-ray. Within diskir, Comptonization is implemented 
with the nth comp model (iZdziarski. Johnson & Magdziarzl 
1 996tZvcki. Done & Smith! 19991) . The physical scenario be¬ 
ing considered is a near-IR to UV thermal component from a 
truncated optically thick accretion disk, providing seed pho¬ 
tons to a Comptonization region with hot electrons. 

The diskir model alone provides a reasonably good de¬ 
scription of the spectrum, but it is not formally acceptable 
with Xv = 1-38 for 2152 dof (see Table0. The fact that the 
thermal component acts as the seed photon distribution for the 
Comptonization emission leads to a somewhat higher value of 
kTm (12!* eV compared to 5! 2 eV found in Section 3.1.2) and 
a lower normalization, corresponding to a somewhat smaller 
disk inner radius. The temperature of the Comptonizing elec¬ 
trons is constrained to be >60keV, and the Comptonizing 
fraction is L c /L c i = 4.2!j {. 

Adding a second continuum component provides a much 
improved fit, and approximately the same improvement is 
seen whether we add a power-law with an exponential cut¬ 
off or a reflection component (see Table 0. Also, both two- 
component models lead to very similar values for the thermal 
component with kT{ n increasing to 29!) 7 eV in one case and 
29 ± 5 eV in the other. The values of Adiskbb decrease further, 
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Fig. 6.— (a) Fit to the Swift J1753.5-0127 spectrum, including the 
data from the near-IR to the X-ray. The model is diskir plus a 
power-law with a high-energy cutoff. The power-law is the harder of 
the two components, (b) Data-to-model ratio. 

but they still imply a large disk truncation radius. 

The two-component models have very different implica¬ 
tions for the properties of the Comptonization region. The 
physical scenario we are considering in adding an extra 
power-law is that either this emission comes from the jet or 
that there is an inhomogeneous or multi-phase Comp t oniza¬ 
tion region (Makis hlma et al.l l2008t iTakahashi et al] 120081 : 
lYamadaetal. I2013fk When this component is added, as 
shown in Figure [6] its best fit parameters imply a very hard 
spectrum T = 1 . 33 ^ 25 , and it dominates at high energies, 
so that the diskir Comptonization component can have 
much lower values of kT e (the constraint is >35 keV) and 
L c /Lj = 0.77 ±0.17. A value of L c /Lj below 1.0 is unusual 
for the hard state, but this is due to the fact that much of the 
hard X-ray flux is in the power-law component. 

On the other hand, when reflection is added, as shown in 
Figure [7] the physical scenario is that the diskir Comp¬ 
tonization component is being reflected from the truncated 
disk. As the diskir component must produce the high en¬ 
ergy emission in this case, a very high Comptonization tem¬ 
perature is required ( kT e > 429 keV) and the Comptonizing 
fraction increases to L c /Lj = 2.4±0.6. The ref lionx pa¬ 
rameters are similar to those described above for the X-ray 
only fits. The ionization state is low, with a value of £ = 5.0 'j'i 
ergcms -1 . Also, the Fe abundance is 0.33 ±0.09, and the 
covering fraction is £l/2n = 0.20. The ionization parame¬ 
ter and the covering fraction do not seem unreasonable for a 
cool and truncated disk, but we cannot say with any certainty 
which two-component model is more likely to be correct. 

3.1.6. Full SED 

When the full radio to hard X-ray SED is put together, it 
is immediately clear that the extrapolation of the power-law 
seen in the radio band is well above the flux measured in the 
near-IR (even after dereddening). This implies that the radio 
component, which is attributed to the compact jet, must have 
a spectral break between the IR and radio bands. Thus, we 
fit the SED with a model consisting of a broken power-law 
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Fig. 7.— (a) Fit to the Swift J1753.5-0127 spectrum, including the 
data from the near-IR to the X-ray. The model is diskir plus a 
reflection component, (b) Data-to-model ratio. 

(bknpower) and a diskir component. The bknpower 

component provides all of the emission in the radio band, and 
we fix the power-law index below the break energy (£b rea k) 
to Ti = 0.7. The index above the break f^) is left as a free 
parameter, and we find that the best fit model has a strong con¬ 
tribution from the bknpower above f«20keV. As described 
above, the GSO data require a cutoff, and we multiplied the 
broken power-law component with a high-energy cutoff. 

The continuum components are multiplied by redden and 
tbabs as described above. We fixed E(B-V) to 0.45 and Ah 
was left as a free parameter. The fit parameters are given in 
Table [6] and the quality of the fit is \ 2 = 1.28 for 2156 dof. 
We left the normalizations between the X-ray instruments as 
free parameters, but we fixed all of the non-X-ray instrument 
normalizations to the SwiftfXRT normalization. 

We used the XSPEC routine steppar to determine the 
range of possible values for Ebreak- The x 2 values are 
nearly constant over a large range, increasing sharply at 
1.0 x 10~ 7 keV (2.4 x 10 10 Hz), which corresponds to the 
highest radio frequency measured, and at 1.5 x lCT 5 keV 
(3.6 x 10 12 Hz). At the upper limit, T 2 becomes steeper to 
avoid over-producing in the near-IR, but \ 2 becomes worse 
because the component no longer extends to the X-ray band. 
For Figure[8j showing the fitted SED, we set Eb rea k to 1.0 x 
10 _6 keV as an example. The main result is that it is pos¬ 
sible for the broken power-law to account for the hard X- 
ray excess. The diskir parameters for the full SED fit 
(see Table [6} are almost the same as the parameters for the 
diskir+highecut*pegpwrlw fits to the near-IR to X- 
ray fits (see Tabled}. 

If the hard X-ray excess is explained by a reflection com¬ 
ponent instead of the broken power-law (i.e., the jet), then Id 
could be steeper, allowing for even higher values of Eb rea k- We 
explored this possibility by fitting just the radio to UV spec¬ 
trum with a modified version of the model shown in Table [6] 
The modifications include removing highecut and fixing 
the diskir components related to Comptonization to the 
values found for the full SED. In addition, while we allowed 
T 2 to be a free parameter, we did not allow this part of the 
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Fig. 8. — (a) The radio to hard X-ray spectral energy distribution for 
Swift J1753.5-0127. The model is a broken power-law with a high- 
energy cutoff (dash-dotted line) and a diskir component, which 
we have divided into its thermal component (dashed line) and its 
Comptonization component (dotted line). The points are not dered- 
dened, and we use E(B-V) = 0.45 and Nu = 2.84 x 10 21 cm' 2 for the 
model. The solid line is the sum of the components, (b) The same 
data and model after dereddening. (c) The same data and model mul¬ 
tiplied by energy. 

broken power-law to be steeper than = 2. While the lower 
limit on £j,reak is unchanged, the upper limit moves higher, 
and values as high as Ebreak = 6.5 x 10 -5 keV (1.6 x 10 l 3 Hz) 
are possible. 


3.2. X-ray Timing 

We made power spectra using the Suzaku/X IS and 
Swift/X RT data. XIS has a larger effective area and the expo¬ 
sure time is much longer than XRT, so the statistical quality is 
much better. However, the XRT data are useful because of the 
higher time resolution. There is good agreement between the 
two power spectra in the frequency region where they over¬ 
lap (see Figure[9]), but precise agreement is not expected due 
to the different times being covered and the slightly different 
energy bandpasses. Thus, we fitted the power spectra sepa¬ 
rately. For XIS, we used a zero-centered Lorentzian and a 
power-law at low frequencies. For XRT, the zero-centered 
Lorentzian is sufficient. The parameters are shown in Table[7] 
and the fractional rms of the Lorentzians are 27.3%±0.2% for 
XIS and 22%±2% for XRT, which is consistent with the rel¬ 
atively high levels of variability expected for the hard state. 
The full-width at half-maximum (FWHM) of the Lorentzians 
are 0.220 ± 0.005 Hz for XIS and 0.33 ±0.07 Hz for XRT. 
In previous work on timing analysis of Swift J1753.5-0127 
dSoleri et alJl201 3b iKalamkar et al.ll2015l) . the Lorentzian fits 
were characterized by the frequency where the power spec¬ 
trum is maximal when plotted as frequency times rms power 
(// m;lx ) as shown in Figure [9] For XIS and XRT, the values of 
t/ max are 0.110±0.003Hz and 0.16±0.04Hz, respectively. 

4. DISCUSSION 



0.001 0.01 0.1 1 

Frequency (Hz) 


Fig. 9.— Soft X-ray power spectrum from Suzaku/X IS (black) and 
Swift/XKT (blue) fitted with a zero-centered Lorentzian and a low- 
frequency power-law. 

In this work, we have performed detailed spectral fits to the 
most complete SED that has been obtained for Swift J 1753.5— 
0127 to date. While previous multi-wavelength studies of this 
source that included radio measurements_have covered the 
radio, near-IR, op t ical, and X-ray dCadpllgJjeLgtaL 2007; 
Durant et all 120091 iRevnolds et al.l l2010fc ISoleri et al. 201C : 
Zhan eWuan & Chatvl 120 1 oFT T we have obtained radio detec¬ 
tions at nine frequencies, included UV coverage, and used 
a combination of seven X-ray spectra, covering 0.5 keV to 
240 keV. Here, we discuss three main topics: 1. the implica¬ 
tions of the constraint on inbreak f° r the compact jet properties; 
2 . what we can infer about the properties of the optically thick 
accretion disk; and 3. the possible origins of the high-energy 
emission components. 

For all these topics, it is useful to estimate the luminosity of 
Swift J1753.5-0127 during these observations. For the model 
shown in Figure[ 8 ] the absorbed flux over the full energy band 
covered (2 x 10 -8 keV to 240 keV) is 1.25 x 10 ' 9 erg cm -2 s' 1 . 
Although there is uncertainty about the break frequency of 
the broken power-law, this leads to very little uncertainty in 
the flux since essentially all of the flux is above 1 eV. The 
unabsorbed flux is 2.71 x 10 _ 9 ergcm _ 2 s _1 in the 1 eV to 
240 keV band, and this represents the bolometric flux. This 
is for the model in Table | 6 j but the unabsorbed flux for the 
diskir+ref lionx model shown in Table [5] gives an un¬ 
absorbed flux of 2.38 x 10 -9 erg cm -2 s -1 due to the lower col¬ 
umn density. Using the average of these two unabsorbed 
fluxes, the bolometric luminosity is 2.7 x 10 36 c/|ergs _1 , 
where d 3 is the distance to the source in units of 3 kpc. For 
a black hole mass of 5 Mr.,, this corresponds to an Eddington- 
scaled luminosity of 0.41% c /3 MJ 1 , where M 5 is the black 
hole mass in units of 5 M 0 . 

4.1. The Compact Jet and the Break Frequency 

We are able to obtain a constraint on //break because of the 
rising and well constrained radio spectrum (a = 0.29 ± 0.05) 
along with the fact that the spectrum rises from A)-band to 
higher frequencies. Without considering the X-rays, we find 
that //break < 1-6 X 10 13 Hz (log //break, Hz < 13-2). If the jet 
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does contribute to the X-rays, then z/b rea k < 3.6 x 10 12 Hz 
(loginbreak,Hz < 12.6). A study of 16 i^reak measurements or 
limits for nine black hole systems in the hard state found 
mostly higher values than the Swift J1753.5-0127 upper lim¬ 
its (iRussell et all I2013al) . For the measurements, the me¬ 
dian value of logi'bieakHz i s 13.68, and the values range 
from 12.65 (for XTE J1118+480) to 14.26 (for GX 339- 
4 and V404 Cyg). When limits are also considered, there 
are still only two measurements that are as low as the value 
found for Swift J1753.5-0127: log ^, re ak Hz = 12.65 ±0.08 for 
XTE J1118+480 and < 13.13 for GX 339-4. 

While relatively low, the single //break measurement for 
Swift J1753.5-0127 does not necessarily indicate any¬ 
thing unusual about the system itself. Multiple mea¬ 
surements of individual systems show significant changes 
for GX 339^1, XTE J1118+480, MAXI J1836-194, and 
MAXI J1659-15 2 (Gand hi et al.ll20TTt IRussell et all20l3allbt 

Ivan der Horst et alJ 120131) ! For GX 339-4. IGandhi et all 

(120111) found that changed by a factor of >10 in less 
than a day. For MAXI J1836—194, six measurements over a 
period of less than two months showed changes in log z/break,Hz 
from close to 11 to close to 14 while the source change d X- 
ray luminosity and hardness (IRussell et al .11201 JR 120141) . and 
the highest value of z/^eak occurred when the source was at its 
lowest X-ray luminosity with its hardest X-ray spectrum. The 
Swift J1753.5-0127 measurements occurred when the spec¬ 
trum was hard and the X-ray luminosity was low; thus, it may 
not follow the same trend as MAXI J1836—194. However, 
this is not surprising since the larger source sample studied in 
IRussell et all (12013al) did not show any evidence for a corre¬ 
lation between X-ray luminosity and r^ rea k- 

In the canonical model for compact jets 
dBlandford & KonigH 119791) . the jet spectrum is com¬ 
posed of a superposition of synchrotron components with 
a continuum of peak frequencies due to changing optical 
depth. The synchrotron spectrum from each region depends 
primarily on the magnetic field strength and also on the 
radial size of the jet. The value of //b rea k depends on both the 
magnetic field and the radial size of the jet in its acceleration 
zone, which is close to the base of the jet. To place constraints 
on these quantities (B and R), we use equations 1 and 2 from 
iGandhi et al.l (1201 It) , which are based on a si ngle-zone 
cylindrical approximation (IChatv. D ubus & Raic hoorll201 ll) . 
We estimate the upper limit on B using the parameters from 
the full SED fit (see Table [Gj. The input parameters to the 
equations are z/^eak < 3.6 x 10 12 Hz, the flux at 3.6 x 10 12 Hz, 
which is 1.42 mJy, and the slope of the power-law above 
z/break■ To determine the slope, we fixed z/b, e ak to 3.6 x 10 12 Hz, 
refit the SED, and found a value of 1.4, which corresponds 
to a = -0.4. The upper limit on the magnetic field strength 
in the acceleration zone is B < 2.4 x 10 3 r/J 024 G and 
R > 1.8 x 10 9 of 9 - 936 cm. If we do not consider the X- 
rays, z^, rea k < 1.6 x 10 I3 Hz, the flux at 1.6 x 10 I3 Hz is 
2.18 mJy, and the slope of the power-law above z^, rea k is 
assumed to be 2 (a = -1), giving B < 9.6 x 10 3 elf 2 ' G and 
R > 4.6 x 10 8 t/ 9 ' 954 cm. Also, from the radio alone, we know 
that z/b r eak > 2.5 x 10 u) Hz, and the flux at this frequency is 
0.34 mJy. Assuming a = -1, we derive B > 18 dj 021 G and 
R< 1.2 x 10 11 4 954 cm. 

Two examples of hard state black hole systems for which B 
and R have been previously calculated using this same tech¬ 
nique are GX 339-4 (IGandhi et al.il20TTl) and MAXI J1836- 


194 (IRussell et all 120141) . For GX 339-4, these quantities 
were estimated to be B ss 1.5 x 10 4 G and R ss 2.5 x 10 9 cm. 
For MAXI J1836—194, estimates for B and R were obtained 
for three hard state observations: one during the rise of an out¬ 
burst and two during outburst decay. Figure 6 of Russ ell et all 
(120141) shows B ~ 10 2 G and R ~ 10 12 cm during the rise and 
B ~ 3 x 10 3-4 G and R ~ 10 9 ~ 10 cm during the decay. Thus, the 
ranges of B= 1.8 x 10‘-9.6 x 10 3 df 21 G and R = 4.6 x 10 8 - 
1.2 x 10 11 a? 9,954 cm that we derive for Swift J1753.5-0127 
are largely consistent with the range of values previously de¬ 
termined for these two sources. The best agreement in the jet 
properties between Swift J1753.5-0127 and GX 339—4 and 
MAXI J1836—194 (during decay) occurs if the actual value of 
z/break for Swift J 1753.5—0127 is close to the upper end of the 
range of possible values. 

4.2. The Optically Thick Accretion Disk 

Here, we discuss the spectral components that can be mod¬ 
eled as thermal emission and the implications for the optically 
thick accretion disk. First, we discuss the near-IR to UV com¬ 
ponent that is consistent with a multi-temperature disk model 
with kT[ n = 28 4 j[ eV. Then, we consider the possibility of a 
second thermal component in the soft X-ray band with k7j n ss 
150 eV. 

Our spectral model assumes that the near-IR to UV emis¬ 
sion is strongly dominated by a disk component, and it is 
worthwhile to consider how secure this assumption is. We 
know that at least a large fraction of the emission comes from 
the disk because of the double-peaked emission lines that 
are seen in this bandpass (Froning et al. 2014; Neustroev 
et al. _2014; Rahoui et al., submitted to ApJ). However, 
iNeustroev et al.l (120141) also find a weak emission line and 
two weak absorption lines (all three unidentified) in the opti¬ 
cal, which they interpret as coming from the companion star. 
If this interpretation is correct (and we note that the fiducial 
black hole mass and source distance that we use in this pa¬ 
per depend on it), then it requires some contribution from the 
companion in the optical. Without X-ray irradiation, the emis¬ 
sion from the companion would be negligible: a b lackbody 
with a temperature of 3000 K ( INeustroev et al.l2014l) . a radius 
equal to the companion’s Roche lobe size of 1.68 x 10 10 cm, 
and a distance of 3 kpc has a flux that is two orders of mag¬ 
nitude lower than the measured flux in the near-IR and three 
orders of magnitude lower in the optical. Thus, the tempera¬ 
ture of the irradiated side of the companion must be signifi¬ 
cantly hotter for there to be a contribution to the optical flux. 
However, the crucial point is that even if the three lines are 
from the companion, they are extremely weak in comparison 
to the very strong double-peaked lines from the disk, indicat¬ 
ing that the disk emission is much stronger than any potential 
contribution from the companion. 

One possibility that we cannot completely rule out is that 
there are additional components from the compact jet. The 
broken power-law emission represents the post-shock syn¬ 
chrotron component. While this is the only component that 
has been seen in SEDs of accreting black holes that is widely 
accepted as emission from the compact jet, theoretical jet 
models indicate that pre-shock synchrotron can be relatively 
bright in the optical and UV dMarkof f. Nowak & Wilrn sl2005l: 
iHoman et al.l 120 05: Migl iari et al.l 12007b iMaitra et al.l 120091) . 
Another possibility that has been suggested as a contributor 
to the optical emission is synchrotron emission from non- 
thermal electrons in the hot accretion flow (i.e., the corona). 
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A complex optical/X-ray cross-co rrelation function was re - 
ported for Swift J1753.5-0127 (iDurant et al.l 120081 1201 ll) . 
and it was shown that it could be explained if the opti¬ 
cal emission had components from the disk and the corona 
dVeledina. Poutanen & Vurml201 ll) . The coronal contribution 
to the cross-correlation function has been observed to vary in¬ 
versely with the strength of the disk (iHvnes et al.ll2009l) . Ra- 
houi et al. (submitted to ApJ) show that Swift J1753.5-0127 
had a strong and likely dominant thermal disk component in 
observations taken a few months after oursFT. which would 
suggest a relatively weak coronal contribution to the optical 
during our observations. 

With the caveats about the possibility of a fractional con¬ 
tribution from pre-shock synchrotron emission or the corona, 
we can compare the parameters of the thermal near-IR to UV 
component to previous studies of Swift J1753.5-0127 SEDs 
wher e this component has also been modeled as therma l emis¬ 
sion (IZhang. Yuan & Chatvll201(llFroning et al.ll2014l) . From 
Tables [3 (last two columns) and [6] the values of /V ( |j.skbb are 
(5.5!^) x 10 7 , (8.2^) x 10 7 , and (9.0+f 9 2 ) x 10 7 . Using 
Equation 1, these values imply strongly truncated disks. As 
before, we assume Mbh = 5M 0 , d = 3 kpc, and i = 40°. To 
determine the lower limits on the inner disk radii, we as¬ 
sume / = 1, and the values are R; n > 227 /A, > 212 R g , and 
> 409 R g . While previous studies have mostly assumed a 
larger distance to Swift J1753.5-0127, this would make the 
values of R m larger. iFroning et al.l (12014l) modeled a near-IR 
to UV SED and determined that R m needed to be > 100 R g 
to avoid overpredicting the simultaneously measured X-ray 
spectrum. iZhang, Yuan & Chatvl (!201 Oi) used a self-consistent 
model with optically thick disk emission, jet emission, and a 
Comptonization component, and they were able to fit_a radio 
to hard X-ray SED with R m = 500 R g . IZhang, Yuan & Chatvl 
(120101) assumed different values for d, Mbh, and i, and if we 
recalculate their R m using the values we adopt, the result is 
Ri n = 350 R g . While the precise value of Rj n is likely to vary 
in time, all of these measurements suggest that the near-IR to 
UV component comes from a strongly truncated disk. 

The spectral fits also constrain the outer disk radius based 
on the parameter log(R out /Ri n ) = 2.3Ujjo4 ( see Tabled- For 
■Vdiskbb = 9 x 10 7 (the best fit value), we calculate R out = 6.6 x 
10 10 cm. We compare this value to the system parameters re¬ 
ported by iNeustroev et all (120141) . where they determine that 
the binary separation is a < 1.1 xlO 11 cm, and the size of 
the black hole’s Roche lobe is 7.1 x 10 I0 cm. A filling frac¬ 
tion of 90% is typically assumed for an accretion disk, which 
would result in a predicted disk size of 6.4 x 10 10 cm, which 
is in excellent agreement with our measurement. Although it 
will be important to confirm the system parameters with ra¬ 
dial velocity measurements of the companion star when the 
source is in quiescence (if it is bright enough), we see this 
R out comparison as another piece of evidence that the near-IR 
to UV component is strongly dominated by emission from the 
accretion disk. 

The 150 eV component is marginally significant in the XRT 
spectrum, and it is not detected when the XIS data are in¬ 
cluded. However, for previous observations of Swift J 1753.5— 
0127, the presence of a 0.1-0.4 keV thermal component 
was well-established from spectral (see references in Section 

25 The Rahoui et al. (submitted to ApJ) observations were made on 2014 
August 16 (MJD 56885). and the X-ray light curves shown in Figure[T|do not 
show any major change between April and August. 


3.1.3) and timing (lUttlev et all 1201 ll) measurements. Even 
though our 2014 April observation is at a moderately lower 
X-ray flux level (only a factor of 2-3 lower than the ma¬ 
jority of the previous observations), seeing a weak ther¬ 
mal component in the X-ray band is not surprising. If we 
use Vdiskbb = 1-1 x 10 4 and carry out the same inner ra¬ 
dius calculation as performed for the near-IR to UV com¬ 
ponent, we obtain R\„ = 5 R g for / = 1 and R m = 14 R„ for 
/ = 1.7, suggesting that this component could come from a 
disk that extends close to the ISCO. The presence of two 
thermal components in the SED of Swift J1753.5-0127 has 
been previously reported dChiang et al.l 12010I) . and poten¬ 
tial physical interpretations are discussed in that work. It 
has been shown that a small inner optically thick accre¬ 
tion disk can form due to condensation of material from 
the corona dLiu et al .1120071 iMever. L iu & Me ver-Hofmeisterl 
]2007l: Taam et al.l 120081) . and Ichiang et al.l (120101) consider 
this possibility as well as a scenario where strong irradiation at 
the inner edge of a truncated disk distorts the temperature pro¬ 
file. For the inner disk possibility, it has been predicted that 
the inner disk can exist down to L/L pm ~ 0. 1% and then com¬ 
pletely evaporate below this level dTaam et al.ll2008l) . Thus, 
given the luminosity of Swift J1753.5-0127 during our ob¬ 
servation (L/L-e&a ~ 0.4%), the presence of an inner disk is 
predicted. 

As previously mentioned, the luminosity at the time of our 
observation of Swift J1753.5-0127 was close to the lowest 
level since the source was discovered, but it was only a fac¬ 
tor of a few times lower than the highest levels seen over the 
past several years (see Figure Q}. The X-ray power spectrum 
also suggests that the properties during our observation were 
at one end of a continuum as opposed to requiring some ma¬ 
jor overall change in the system. ISoleri et akl (i2013l) report 
on timing analysis of 67 Rossi X-ray Timing Explorer (RXTE) 
observations of Swift J1753.5-0127 during 2009 and 2010. 
While the comparison to our observation is somewhat com¬ 
plicated by the fact that most of the RXTE observations re¬ 
quired two Lorentzian components, fifteen of the power spec¬ 
tra were fitted with a single Lorentzian, allowing for a direct 
comparison. For those cases, the values of // max range from 
0.18 Hz to 3.18 Hz. Thus, our Suzaku and Swift measurements 
of 0.110 ±0.003 Hz and 0.16 ±0.04Hz, respectively, are only 
slightly lower than the lSoleri et al .l ( 201 31) measurements. 

4.3. The Origin of the X-ray Emission 

A major question in recent years concerns how much of 
the X-ray emission can be attributed to the compact jet. 
In the model of iMarkoff. Nowak & Wilmsl (120051) . the jet 
can produce X-rays via post-shock synchrotron emission, 
which can be modeled as the broken power-law that we 
use in our fits, or synchrotron self-Compton (SSC) from 
the base of the jet, which can contribute in the hard X- 
ray band. The SEDs of GX 339-4, GRO J1655-40, and 
XTE J1118+480 allow for the possibility that all the soft X- 
ray em ission comes from the post-shock synchrotron com¬ 
ponent (IM arkoff. Nowak & Wilmsl2 005: Mi gliari et alJ2007t 
iMaitra et al .1120091) . For Swift J1753.5-0127. Figure [8] shows 
that such a scenario is ruled out, and a Comptonization com¬ 
ponent is strongly required by the data. 

The question of what makes Swift J1753.5-0127 different 
is directly relevant to the question of what is different about 
the outliers in the X-ray/radio correlation. Although one 
possibility is that Swift J1753.5-0127 has a stronger Comp¬ 
tonization component in the X-rays, another possibility is that 
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it has a weaker radio jet. In Section 4.1, we showed that the 
highest possible peak flux for the Swift J1753.5-0127 bro¬ 
ken power-law component is 2.18 mJy. This corresponds to 
a specific (peak) luminosity of 2.3 x 10 19 t / 2 erg s“'Hz~' at 
J'break = 1 .6 x 10 13 Hz. iRussell et al.i ( 2013a ) give 15 peak 
luminosities for nine hard state black hole systems, and the 
values range from 7.1 x 10 19 to 1.9 x 10 22 erg s _l Hz -1 with a 
median value of 1.2 x 10 21 erg s -1 Hz -1 . Thus, assuming a dis¬ 
tance of 3 kpc, the peak jet luminosity for Swift J1753.5-0127 
is 50 times lower than the median and 3 times lower than the 
least luminous system. The distance to Swift J1753.5-0127 
would need to be 5-6 kpc to to move the Swift J1753.5-0127 
peak jet luminosity close to the least luminous system, which 
is conceivable, but it would need to be r ;21 kpc to make the 
Swift J1753.5-0127 comparable with the median peak jet lu¬ 
minosity, which can be ruled out. 

We made a second radio luminosity comparison by inte¬ 
grating the radio power-law measurements for Swift J 1753.5— 
0127 and the black hole sources from IRussell et al] (12013al) 
up to zzbreak- For Swift J1753.5-0127, the luminosity up 
to 1 .6 2<. 10 i3 Hz is 3.1 x 10 32 c/f ergs~'. For the sources 
from IRussell et alj (1 2013ai) . not all 15 of the SEDs are high 
enough quality to make a reliable luminosity determination. 
There was sufficient information to calculate ten luminosi¬ 
ties for eight sources. These ranged from 1.8 x 10 33 ergs _1 
for XTE J1118+480 and 2.0 x lO^ergs " 1 for Cyg X-l to 
1.1 x 10 36 ergs -1 for GS 1354-64 and 3.1 x lO^ergs ” 1 for 
V404 Cyg. The median value is 1.1 x 10 35 ergs _1 , and the 
distance to Swift J1753.5-0127 is certainly not large enough 
for the luminosity to approach that value. Thus, the low lu¬ 
minosity radio jet may be at least part of the reason why 
Swift J1753.5-0127 is an outlier. 

While the Swift J1753.5-0127 SED is consistent with 
Comptonization being dominant at soft X-rays, our results 
show that multiple components are required to explain the 
entire 0.5-240 keV X-ray spectrum. In our spectral fits, we 
considered a reflection component or the post-shock syn¬ 
chrotron component. Figure [7] illustrates the reflection pos¬ 
sibility, and such a scenario is consistent with our overall 
picture for the system. The outer optically thick disk could 
produce a weak (fl/27r = 0.20) reflection component, and it 
would be expected to have a low ionization, which is consis¬ 
tent with £ = 5.0^ 2 erg cm s _I . While an iron line detection 
would be strong evidence in favor of the reflection interpreta¬ 
tion, there is no iron line in the Swift J1753.5-0127 spectrum, 
but we find that it is possible to explain the lack of an iron 
line with an iron abundance of 0.33 ± 0.09 of the solar value. 
This iron abundance may be problematic for the reflection in¬ 
terpretation, but we do not think that it is low enough to rule 
it out. We have also considered the fact that this is the only 
model that requires a very high Comptonization temperature 
(kT e > 429 keV). This occurs because the reflection compo¬ 
nent falls at high energies, allowing the overall model to fit the 
steeper Suzakul GSO spectrum without an exponential cutoff 
in the direct model. This electron temperature is higher than 
has been inferred from measurements of other accreting black 
holes in the higher luminosity parts of their hard states, which 
are typically in the 50-120 keV range (iPoutanen & Veledinal 
120141 . and references therein). However, it is predicted that 
kT e should increase to hun dreds of keV in the low er luminos¬ 
ity parts of the hard state dGardner & Donell2013l) . Thus, the 
lack of an iron line is the strongest reason to disfavor the re¬ 
flection possibility, but this scenario is not ruled out. 


The model where the hard X-rays are due to post-shock 
synchrotron emission (see Figure [ 6 ]> has the advantage of a 
much more typical electron temperature (kT e > 33 keV). On 
the other hand, the slope of the power-law, T = 1 ,337q 25 
(a = -0.33^'gg), while not unreasonable for optically thin 
synchrotron emission, is harder than is seen for other black 
hole systems, which have values of a between - 0.68 and - 
1.38 (IRussell et al.ll2013al) . Such a hard spectrum also re¬ 
quires that the spectrum is sharply cutoff to explain the steeper 
Suzakul GSO spectrum, and the exponential cutoff with E cul = 
20 ± 3 keV and £f 0 id = 142£'g° keV is probably inconsistent 
with the more gradual cutoff predicted for a synchrotron spec¬ 
trum (IZdziarski et al] 120031) . A third possibility that was 
mentioned above but not specifically considered in our spec¬ 
tral modeling is that the extra hard X-ray component is due 
to SSC emission from the _base of the jet. Fits with the 
iMarkoff, Nowak & Wilmsl ( 2005 ) compact jet model are be¬ 
yond the scope of this paper, but it would be interesting to use 
our SED to test this model in future work. 

Finally, we have considered whether any of our conclusions 
might be affected by day-to-day source variability given that 
the observations we use for the full SED cover ~2.8 days 
from the X s -band observation to the VLA observation (al¬ 
though most of the measurements for the SED come from 
a smaller span of times). Figure U shows that there is little 
day-to-day variability in the optical and near-IR during the 
campaign, and this is consistent with previous long-term stud¬ 
ies of Swift J1753.5-0127 (e.g JShaw et aLi Z OO) . Also, the 
Suzaku/X IS observations show day-to-day stability in the soft 
X-ray flux. It is a little less clear whether there are changes in 
the radio and IR compact jet spectrum as other black hole sys¬ 
tems have shown significant changes in the break frequency 
on time scales of a day as discussed in Section 4.1. We al¬ 
ready consider a large range of break frequencies; thus, the 
conclusions that there is a separate thermal component in the 
near-IR to UV and that a Comptonization component is re¬ 
quired in the soft X-ray should not be affected. However, the 
question of whether the extra hard X-ray component comes 
from the compact jet depends very sensitively on the break 
frequency and spectral slope. To reach a definitive conclusion 
on the origin of the hard X-ray emission may require simulta¬ 
neous radio and hard X-ray monitoring. 

5. SUMMARY AND CONCLUSIONS 

We have obtained radio, near-IR, optical, UV, and 
X-ray coverage for the long-term black hole transient 
Swift J1753.5-0127 in 2014 April when the source was in 
the hard state at one of its lowest X-ray luminosities (2.7 x 
10 36 d\ erg s ~ 1 j since the discovery of the source. We per¬ 
formed fitting of the broadband energy spectrum as well as 
the X-ray power spectrum. We obtain results concerning the 
compact jet, the optically thick accretion disk, and the origin 
of the X-ray emission, which is also relevant for the question 
of why Swift J1753.5-0127 is a radio/X-ray correlation out¬ 
lier. 

With the combination of the rising radio spectrum, and the 
rise in the near-IR, r^ rea k is constrained for the post-shock syn¬ 
chrotron component of the compact jet, and this provides con¬ 
straints on B and R for the jet acceleration zone. While the 
post-shock synchrotron component may contribute in hard X- 
rays, the soft X-ray flux is far too high to be part of this com¬ 
ponent, which we model with a Comptonization component. 
Based on this result. Swift J1753.5-0127 appears to be an out- 
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lier because of the combination of a strong Comptonization 
component and a jet with peak and broadband luminosities 
significantly lower than is seen for other black hole systems. 

The low jet luminosity and the low extinction for 
Swift J 1753.5-0127 appear to provide an opportunity to 
clearly see emission components that may be too weak or too 
absorbed to see in other systems. The double-peaked emis¬ 
sion lines (Froning et al. 2014; Neustroev et al. 2014; Rahoui 
et al., submitted to ApJ) clearly show that the near-IR to UV 
spectrum has at least a strong (likely dominant) thermal disk 
component. Further evidence that the near-IR is dominated 
by thermal disk emission is that the component can be mod¬ 
eled by a disk with an outer radius of R ou t/Rg = 90,000 d /3 MJ 1 
(R out = 6.6 x 10 10 d /3 cm), consistent with the expected size of 
the disk given previous measurements of the size of the com¬ 
panion’s Roche lobe. The fact that this component does not 
contribute in the X-ray band constrains the inner radius to be 
Rin/Rg > 212 d /3 M 5 '. While this implies that the near-IR to 
UV emission comes from a strongly truncated disk, there is 
also some evidence for a weak 150 eV thermal component in 
the soft X-rays, and its inner radius could be as small as 5- 
14 R g . The presence of two thermal components could pro¬ 
vide support for predictions that low luminosity systems may 
have inner and outer optically thick disks with a gap in the 
middle. 

Finally, we have considered the possibility that there is a 
reflection component in the spectrum. In the presence of 
strong hard X-rays, one expects to see a reflection compo¬ 
nent from the optically thick material. The hard X-ray spec¬ 
trum is consistent with the presence of a reflection compo¬ 
nent, but no iron line is detected. The low ionization (£ = 
5.0^4 erg cm -2 s -1 ) and low covering fraction {VL/2-k = 0 . 2 ) 
would favor the possibility that this component comes from 
the outer optically thick disk. If reflection is the cause of the 
second hard X-ray component, then invoking the jet to ex¬ 
plain the extra hard X-ray emission (see Figure [ h} may not be 
required. 
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TABLE 1 

Observing Log and Exposure Times 


Mission 

Instrument 

Energy/Filter 

ObsID 

Start Time (UT) 
(in 2014) 

End Time (UT) 
(in 2014) 

Exposure 

(S) 

Radio 

VLA 

— 

4-8 GHz 

— 

Apr 5, 11.00 h 

Apr 5, 13.00 h 

1518 

VLA 

— 

18-26 GHz 

— 

” 

” 

1746 

AMI 

— 

12-17.9 GHz 

— 

Apr 4, 3.07 h 

Apr 4, 7.56 h 

16,164 

AMI 

— 

12-17.9 GHz 

— 

Apr 5, 2.67 h 

Apr 5, 7.64 h 

17,892 

Near-IR to UV 

Kanata 

HONIR 

BIVIJIKs 

— 

Apr 2, 17.3 h 

Apr 7, 19.2 h 

see text and Table [2] 

SLT 

U42 

g'/r'/i'/z' 

— 

Apr 3, 9.1 h 

Apr 5, 11.0 h 

180^ 

Swift 

UVOT 

V 

00080730001 

Apr 5, 0.4 h 

Apr 5, 5.4 h 

182 

Swift 

UVOT 

b 

00080730001 

” 

” 

182 

Swift 

UVOT 

u 

00080730001 

” 

” 

182 

Swift 

UVOT 

uvw 1 

00080730001 

” 

” 

364 

Swift 

UVOT 

uvm2 

00080730001 

” 

” 

591 

Swift 

UVOT 

uvw2 

00080730001 

” 

” 

731 




X 

ray 



Swift 

XRT 

0.5-1 OkeV 

00080730001 

Apr 5, 0.4 h 

Apr 5, 5.4 h 

2,372 

Suzaku 

XIS0/1/3 

1.2-12 keV 

409051010 

Apr 3, 17.65 h 

Apr 5, 10.69 h 

59,711 

NuSTAR 

FPMA/B 

3-79keV 

80002021003 

Apr 4. 21.35 h 

Apr 5, 12.69 h 

61,038 

Suzaku 

HXD/PIN 

13-65 keV 

409051010 

Apr 3, 17.65 h 

Apr 5, 10.69 h 

50,434 

Suzaku 

HXD/GSO 

50-240 keV 

409051010 

Apr 3, 17.65 h 

Apr 5, 10.69 h 

50,434 


“Each exposure listed in Table[3]was 180 s. 


TABLE 2 

Log of Kanata/HONIR Observations 


Epoch 

Filter 

MJD 

Exposure (s) 

i 

V 

56749.7221 

136 

i 

J 

56749.7232 

120 

i 

K s 

56749.7302 

60 

i 

B 

56749.7405 

75 

3 

J 

56751.7170 

120 

3 

V 

56751.7173 

136 

5 

J 

56754.8052 

120 

5 

V 

56754.8059 

136 
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TABLE 3 

Log of SLT/U42 Observations 



Times of g 

Times of r 

Times of i 

Times of z 


Exposures 

Exposures 

Exposures 

Exposures 

Epoch 

(MJD-56750) 

(MJD-56750) 

(MJD-56750) 

(MJD-56750) 

2 

0.3789 

0.3806 

0.3822 

0.3839 

2 

0.3856 

0.3872 

0.3889 

0.3905 

2 

0.3922 

0.3938 

0.3955 

0.3971 

2 

0.3987 

0.4004 

0.4020 

0.4037 

2 

0.4054 

0.4070 

0.4087 

0.4103 

2 

0.4120 

0.4136 

0.4153 

0.4169 

2 

0.4186 

0.4202 

0.4219 

0.4235 

2 

0.4252 

0.4268 

0.4285 

0.4301 

2 

0.4318 

0.4334 

0.4351 

0.4367 

2 

0.4384 

0.4400 

0.4417 

0.4433 

2 

0.4450 

0.4466 

0.4483 

0.4499 

2 

0.4516 

0.4532 

0.4549 

0.4566 

3 

1.7601 

1.7578 

1.7624 

1.7647 

3 

1.7695 

1.7672 

1.7718 

1.7741 

3 

1.7789 

1.7766 

1.7812 

1.7835 

3 

1.7883 

1.7860 

1.7906 

1.7929 

3 

1.7977 

1.7954 

1.8000 

1.8023 

3 

1.8071 

1.8048 

1.8094 

1.8117 

3 

1.8165 

1.8142 

1.8188 

1.8211 

3 

1.8259 

1.8236 

1.8282 

1.8305 

3 

1.8354 

1.8330 

1.8377 

1.8400 

3 

1.8450 

1.8428 

1.8474 

1.8498 

3 

1.8545 

1.8523 

1.8569 

1.8592 

3 

1.8641 

1.8617 

1.8664 

— 

4 

2.3799 

2.3816 

2.3833 

2.3849 

4 

2.3866 

2.3882 

2.3899 

2.3915 

4 

2.3932 

2.3948 

2.3965 

2.3981 

4 

2.3998 

2.4014 

2.4031 

2.4047 

4 

2.4064 

2.4081 

2.4097 

2.4114 

4 

2.4130 

2.4147 

2.4163 

2.4180 

4 

2.4196 

2.4213 

2.4229 

2.4245 

4 

2.4262 

2.4279 

2.4295 

2.4311 

4 

2.4328 

2.4345 

2.4361 

2.4378 

4 

2.4394 

2.4411 

2.4427 

2.4444 

4 

2.4460 

2.4477 

2.4493 

2.4510 

4 

2.4526 

2.4543 

2.4560 

2.4576 
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TABLE 4 

Parameters for X-ray Spectral Fits 


Parameter 

Units 

pegpwrlw 

highecut* 

pegpwrlw 

highecut* 

pegpwrlw! 

diskbb 

highecut* 
pegpwrlw! 
reflionx 

highecut* 

pegpwrlw! 

diskbb! 

reflionx 

A h 

r 

Flux^ 

£cut 

Cfold 

klin 

^diskbb 

c 

Fe/solar 

^refl 

n/2n 

10 21 cnT 2 
Photon index 

10 12 erg citT 2 s 1 
keV 

keV 

keV 

erg cm s 1 

10 + 

2.01 ±0.05 
1.722 ±0.003 

253 ±3 

2.01 ±0.06 
1.721 ±0.003 

253 ±3 
66 !‘ s 

217 +f 2 

2.35 ±0.11 
1.699 ±0.005 

253 ±4 

6 o;l 2 

209 + g® 8 

0-67+8:81 

3-s+Jf 

2.60 ±0.09 
1.774 ±0.006 

254 ±3 

>411 

<5.3 

0.28 ±0.08 

4 1+0.5 
^* a -3.3 

0.20 

2.51 ±0.13 
1.738 ±0.015 

253 ±4 

67+2° 
o'_i2 
-,61+218 
- 01 -95 

0.67 ±0.06 

2 1+ 1 - 7 

<11 

n 47+0-21 

u -+'-0.15 

9 9+1.0 
-2.0 

0.12 

CxRT 

— 

1.0 

1.0 

1.0 

1.0 

1.0 

CxiS03 

— 

1.06 ±0.01 

1.06±0.01 

1.04 ±0.01 

1.04 ±0.01 

1.04 ±0.01 

Cxisi 

— 

0.97 ±0.01 

0.97 ±0.01 

0.95 ±0.01 

0.95 ±0.01 

0.95 ±0.01 

Cfpma 

— 

1.08 ±0.01 

1.08 ±0.01 

1.07 ±0.01 

1.06 ±0.01 

1.06±0.01 

Cfpmb 

— 

1.10 ±0.01 

1 .10±0.01 

1.08 ±0.01 

1.08±0.01 

1.08 ±0.01 

Cpin 

— 

1.31 ±0.02 

1.31 ±0.02 

1.26 ±0.02 

1.25 ±0.01 

1.25 ±0.02 

Cgso 

— 

1.05 ±0.06 

1.19 ± 0.08 

1.14 ±0.09 

1.16 ±0.07 

1.17 ±0.08 

X 2 /dof 

— 

2990/2143 

2970/2141 

2751/2139 

2715/2138 

2699/2136 


“unabsorbed 2-10 keV, power-law only 


TABLE 5 

Parameters for Near-IR, Optical, UV, and X-ray Spectral Fits 


Parameter 

Units 

diskir 

diskir+ 

diskir+ 




highecut* 

pegpwrlw 

reflionx 

E(B-V) 

— 

0.45^ 

0.45 fl 

0.45 fl 

N h 

10 21 cm -2 

2.08 ±0.05 

2.83 + 0.10 

2.60 + 0.08 

diskir 

t-7in 

eV 

12 +8 

29 + ) 7 

29 ±5 

Ndiskbb 

10 7 

q'7+607 

87 -36 

8-2l 2 , 

5 5 +18 

j.j_ 3 i 

r 

Photon index 

1.734 ±0.003 

1.90 + 0.07 

1.777 + 0.006 

kT c 

keV 

>60 

>35 

>429 

L c /L d 

— 

4 9+2-4 

0.77 + 0.17 

2.4 ±0.6 

/in 

— 

0.1 fl 

0.1° 

0.1° 

Hit 

— 

1.1° 

1.1° 

1.1° 

/out 

— 

0.0“ 

0.0" 

0.0" 

log rout 


i q-j+0.06 

1.85-0 40 

a OO+0.29 
-7—0.11 

9 <0+0-34 

highecut*pegpwrlw 

r 

2nd Photon index 

— 

l oo+o.os 

— 

Flux^ 

1(T 12 erg cm -2 s -1 

— 

<o+35 

6o_42 

— 

£cut 

keV 

— 

2X% 

— 

Cfold 

keV 

— 

151 +63 

— 

reflionx 

? 

erg cm s 1 

— 

— 

5.0/+“ 

Fe/solar 

— 

— 

— 

0.33 + 0.09 

Ffold 

keV 

— 

— 

>507 

/Vrefl 

1CU 1 

— 

— 

0.78/°;l| 

Q,/2tt 

— 

— 

— 

0.20 

CxRT 

— 

1.0 

1.0 

1.0 

CxiS03 

— 

1.06 ±0.02 

1.04 + 0.01 

1.04 + 0.01 

Cxisi 

— 

0.97 ±0.01 

0.94 + 0.01 

0.95 + 0.01 

Cfpma 

— 

1.08 ±0.01 

1.06 + 0.01 

1.07 + 0.01 

Cfpmb 

— 

1.09 ±0.01 

1.07 + 0.01 

1.08 + 0.01 

Cpin 

— 

1.30 ±0.02 

1.25 + 0.02 

1.26 + 0.02 

Cgso 

— 

i i -t+0.06 

1 ‘ 1 '-0.03 

1.16 + 0.07 

1.23 + 0.07 

X 2 /dof 

— 

2963/2152 

2765/2148 

2769/2148 


“Fixed. 

b unabsorbed 2-10 keV, power-law only 
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TABLE 6 

Parameters for SED Fits 


Parameter 

Units 

Value 

E(B-V) 

— 

0.45^ 

Na 

10 21 cm- 2 

2.84 ±0.11 

diskir 

Min 

eV 

28«; 

^diskbb 

10 7 

9.0+f 2 

r 

Photon index 

1.90 ±0.02 

kTe 

keV 

>33 

Lc / L d 

— 

rv q--) +0.04 

U ’ S2 -0.19 

,/in 

— 

0 .1° 

n rr 

— 

i.i a 

/out 

— 

0 .0“ 

log Tout 

log (Rout/Kin) 

9 9 1 +0.06 
-0.04 

highecut*bknpower 

r. 

Below -Ebreak 

0.7“ 

r 2 

Above Ebreak 

i 99+O.IO 

1 • ^ -0.13 

^break 

1(T 6 keV 

0.1-15 

^break 

Hz 

2.4 x 10 10 —3.6 x 10 12 

Normalization 

ph cm -2 s -1 keV -1 at 1 keV 

64 ±2 

Ecut 

keV 

20 ±3 

Efdd 

keV 

142+‘‘° 

CxRT 

— 

1.0 

CXIS03 

— 

1.04 ±0.01 

Cxisi 

— 

0.94 ±0.01 

Cfpma 

— 

1.06 ±0.01 

Cfpmb 

— 

1.07 ±0.01 

Cpin 

— 

1.25 ±0.02 

Cgso 

— 

1.16 ±0.07 

X 2 /dof 

— 

2768/2156 


“Fixed. 


TABLE 7 

Parameters for Power Spectrum Fits 


Parameter 

Units 

Value 

Suzaku (Lorentzian plus power-law) 

^max 

Hz 

0.110 ±0.003 

nns Lor 

— 

27.3%±0.2% 

Power-law index — 

2.12±0.05 

rmSp\ 

0.0001-1 Hz 

4.4%±1.4% 

X 2 /dof 

— 

284/191 

Swift (Lorentzian) 

^max 

Hz 

0.16±0.04 

rms Lor 

— 

22 %±2% 

X 2 /dof 

— 

51/42 





















